Lecture  4:  Local  stretching  theories 


Jean-Luc  Thiffeault 

Notes  by  David  Goluskin  and  Anubhab  Roy 
24  June  2010 


1  Introduction 

In  this  lecture  we  will  try  to  understand  local  or  Lagrangian  theories  involved  in  mixing 
of  passive  scalars.  This  involves  solving  the  advection-diffusion  (AD)  equation  along  fluid 
trajectories.  The  origins  of  these  local  theories  can  be  traced  to  Batchelor’s  idea  of  describing 
the  flow  via  spatially-constant  strain-rate  matrices  with  prescribed  time  dependence  [2]. 
Kraichnan  addressed  the  problem  next  by  considering  the  velocity  field  in  the  AD  equation 
to  be  a  stochastic  Gaussian  field  with  a  time  correlation  that  decays  infinitely  rapidly  (or 
is  white  in  time)  and  a  spatial  correlation  that  has  a  power-law  structure  [6].  This  led 
to  solving  a  stochastic  differential  equation.  Zeldovich  encountered  the  problem  in  the 
context  of  heat  diffusion  and  the  magnetic  dynamo  and  adopted  a  random  matrix  theory 
approach  [13-15].  More  recently,  tools  from  large  deviation  theory  and  path  integration 
have  aided  in  obtaining  a  complete  solution  of  the  problem,  as  will  be  discussed  in  this 
lecture. 

Let  us  revisit  the  AD  equation, 

dte  +  u  ■  VO  =  /tV2(9,  n  =  M2  or  M3.  (1) 

As  mentioned  before  6(x,t)  denotes  a  passive  scalar  field  (say  concentration  of  dye)  and  u(x,t ) 
is  the  ambient  velocity  field  which  stirs  the  scalar  as  it  diffuses  with  a  diffusivity  k. 

Previously  (Lecture  1)  we  have  examined  how  a  “patch”  of  dye  evolves  in  an  extensional 
flow,  u  =  ( Xx ,  —  Ay).  We  will  now  address  the  evolution  of  a  passive  scalar  when  acted  upon 
by  an  arbitrary  linear  flow.  By  definition,  for  a  linear  flow  the  velocity  field  depends  linearly 
on  the  position.  Rigid  body  rotation,  extensional  flow,  and  simple  shear  are  some  classic 
examples  of  linear  flow.  In  general, 

u(x,t)  =  U(t)  +  x  ■  A(t),  Tr(A)  =  0,  (2) 

where  A(t )  is  the  velocity  gradient  tensor.  The  requirement  for  it  to  be  traceless  is  a  conse¬ 
quence  of  incompressibility  of  the  flow  field  (V  •  u  =  0).  For  turbulent  or  Lagrangian-chaotic 
flows,  A(f)  is  a  random  matrix,  having  a  finite  correlation  time,  which  is  the  Lagrangian 
correlation  time  of  the  velocity.  Here  U ( t )  is  a  uniform  flow.  The  advantage  of  using  linear 
flows  lies  in  making  the  problem  analytically  tractable  and  most  importantly  it  serves  as  a 
nice  prototype  for  more  complicated  flows.  We  will  simplify  things  further  by  considering 


34 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

24  JUN  2010 

4.  TITLE  AND  SUBTITLE 

Local  Stretching  Theories 


6.  AUTHOR(S) 


2.  REPORT  TYPE 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Woods  Hole  Oceanographic  Institution, Woods  Hole, MA, 02543 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

00-00-2010  to  00-00-2010 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADA544302 

14.  ABSTRACT 

In  this  lecture  we  will  try  to  understand  local  or  Lagrangian  theories  involved  in  mixing  of  passive  scalars. 
This  involves  solving  the  advection-di  usion  (AD)  equation  along  uid  trajectories.  The  origins  of  these  local 
theories  can  be  traced  to  Batchelor’s  idea  of  describing  the  ow  via  spatially-constant  strain-rate  matrices 
with  prescribed  time  dependence  [2].  Kraichnan  addressed  the  problem  next  by  considering  the  velocity 
eld  in  the  AD  equation  to  be  a  stochastic  Gaussian  eld  with  a  time  correlation  that  decays  in  nitely  rapidly 
(or  is  white  in  time)  and  a  spatial  correlation  that  has  a  power-law  structure  [6].  This  led  to  solving  a 
stochastic  di  erential  equation.  Zeldovich  encountered  the  problem  in  the  context  of  heat  di  usion  and  the 
magnetic  dynamo  and  adopted  a  random  matrix  theory  approach  [13{15].  More  recently,  tools  from  large 
deviation  theory  and  path  integration  have  aided  in  obtaining  a  complete  solution  of  the  problem,  as  will  be 
discussed  in  this  lecture. 


15.  SUBIECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

12 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


(a)  (b) 


Figure  1:  Advection  of  a  concentration  field  by  a  simple  shear  flow,  u  =  (Ux,  Uy )  =  (7?/,  0) 
(neglecting  diffusion,  7  assumed  constant).  An  intially  isotropic  concentration  field  (a)  gets 
both  stretched  and  rotated  by  the  background  flow  (b). 

the  limit  of  high  Schmidt  number,  which  is  a  dimensionless  quantity  defined  as 

Sc  :=  v / k 

where  v  is  the  kinematic  viscosity  of  the  fluid  and  k  as  already  mentioned  is  the  diffusivity 
of  the  scalar.  (If  the  scalar  is  heat,  then  the  Schmidt  number  becomes  the  Prandtl  number.) 
Momentum  diffuses  significantly  faster  than  the  passive  scalar  for  flows  with  large  Sc.  This  is 
the  desired  limit  for  chaotic  advection:  Sc  3>  1  ensures  that  the  velocity  field  appears  much 
smoother  than  the  passive  scalar  concentration,  allowing  the  linear  local  approximation  (2). 
This  is  the  regime  that  was  studied  by  Batchelor  and  leads  to  the  celebrated  Batchelor 
spectrum.  For  Sc  <  1  we  have  the  Kolmogorov-Obukhov-Corrsin  regime,  where  the  scalar 
gets  advected  by  the  turbulent  inertial  range,  unlike  the  Batchelor  regime.  (We  will  not  be 
discussing  the  KOC  regime  in  this  lecture.) 

2  ‘Inertia  tensor’  of  a  patch  of  passive  scalar 

It  is  of  great  interest  to  understand  the  response  of  a  passive  scalar  to  a  flow  with  chaotic 
Lagrangian  trajectories.  As  already  discussed,  a  decent  approximation  to  such  a  velocity 
field  can  be  made  by  superposing  random  linear  flows  (in  the  limit  Sc  >  1).  In  order 
to  get  a  statistical  description  of  the  passive  scalar  held  subjected  to  such  an  ensemble  of 
linear  hows,  we  need  to  understand  the  behaviour  of  the  moments  of  9{x,t).  We  use  angle 
brackets  to  denote  the  volume-integrated  value  of  a  quantity  (similar  to  an  average), 

(/)=  [  fdV. 

J  Q 

We  have  already  seen  (Lecture  1)  that  when  the  domain  is  either  periodic  or  is  bounded  by 
insulating  walls  then 

dt(9)  =  0,  (3) 
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that  is,  the  mean  value  of  9  is  conserved.  This  is  a  consequence  of  the  no-net-flux  boundary 
condition  imposed  on  the  system.  To  track  the  evolution  of  ‘a  blob  of  dye,’  one  can  define 
the  ‘center  of  mass’  of  9  as 

_  (xtf) 

Ci~  (0) 

To  obtain  an  evolution  equation  for  c  we  consider  the  first  moment  of  the  AD  equation, 
and  integrate  by  parts  a  few  times: 


dt(xi0)  +  (XiV  •(([/  +  x  •  A)9))  =  K(XiV29 ), 

=>•  dt(xi9)  -  (( Uj  +  xiAej)9djXi )  =  0, 

=>-  dt{xi9)  -  Ui{9)  -  Agi(d)c£  =  0,  (djXi  =  (%) 


which  leads  to 


c  =  U  +  c  •  A. 


(5) 


Here  the  overdot  represents  a  derivative  with  respect  to  time.  In  the  calculation  the  contri¬ 
bution  due  to  the  dissipative  term  disappeared  due  to  the  no-net-flux  boundary  conditions. 
Thus  the  center  of  mass  of  the  patch  gets  advected  by  the  background  flow. 

Proceeding  to  the  calculation  for  the  next  moment  we  define  the  ‘inertia  tensor’  of  the 
patch  as 


Mij  = 


{ XiXj9 ) 

~~ w  CiCJ ' 


The  additional  —CiCj  term  makes  M  translationally  invariant.  By  construction  M  is  also 
a  symmetric  matrix.  We  proceed  in  similar  fashion  as  above  to  obtain  the  corresponding 
evolution  equation  for  M. 


dt(xiXj9 )  +  (xiXjV  ■  (( U  +  x  ■  A) 9))  =  n{xiXjW29) 

Oj  (xjXj  9  j  Uj(xi9)  Uj  ( x  j  9 )  Apj(xiXp9)  Apt  \  X-jXpO)  —  2KSpj  (9'j 

—t"  (h  T  CiCj )  Ujd  UiCj  -Apj  ( Mip  T  c} cp j  Api  (Adjp  T  CjCp )  —  2 KSij 

=t>  Oj  Mij  -)-  Cj  ( Of  Cj  Uj  CpApj)  -\-  Cj  (Otd  Ui  CpApi}  =  MipApj  T  ApiMjp  T  2  k5^ 


so  that  finally 

M  =  M  ■  A  +  A]  ■  M  +  2kI  .  (6) 

Thus  we  have  obtained  an  evolution  equation  for  the  ‘inertia  tensor’  of  the  patch.  Here  A t 
denotes  the  transpose  of  the  velocity  gradient  tensor. 

The  solution  for  c(t )  and  M  ( t )  can  be  explicitly  written  as 

.  pt 

c(t)  =  c( 0)  •  efo  A(rOdTi  +  /  [/(T)  .  eU  A(n)dndr  ?  (7) 

Jo 

r>{ 

M(t )  =  eJo  Af(n)dn  .  +2f,  gf*  At(n)dn  .  e/*  Mn)dndr  .  (8) 

Jo 

Here  the  quantity  eB  is  a  matrix  exponential,  defined  by  the  usual  power  series 

eB  =  I  +  B +  \B  B  +  .... 
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It  should  be  noted  that  in  the  expression  (8)  for  M(t )  the  two  matrix  exponentials  cannot 
be  combined  and  written  as  efT(A^+A)(Tl')dT1 .  For  matrix  exponentials,  eBl+B2  =  eBl  ■  eB2 
only  if  the  matrices  B\  and  B2  commute.  The  velocity  gradient  matrix,  A,  is  typically  a 
non- normal  matrix  (AA^  7^  A^  A). 

For  a  simple  shear  flow,  we  have  u  =  (jy,  0)  and 


A  = 


0 

7 


Since  A  is  nilpotent  ( A 2 
representation,  so  that 


0),  the  matrix  exponential  has  only  two  terms  in  its  power  series 


eAt  =  I  +  At  = 


1  0 

jt  1 


3  Computation  of  Lyapunov  exponents 

Since  M  is  a  symmetric  matrix  the  eigendecomposition  takes  the  simple  form 

M  =  RDR]  (9) 

where  R  is  an  orthogonal  matrix  (the  columns  of  which  are  eigenvectors  of  A/),  and  D  is 
real  and  diagonal  (having  the  eigenvalues  of  M  on  the  diagonal).  When  the  eigenvalues  are 
sufficiently  distinct,  the  radial  and  angular  degrees  of  freedom  decouple  from  each  other,  as 
we  will  now  show.  (Our  derivation  parallels  that  of  [1].)  Let  us  substitute  the  expression 
(9)  in  the  evolution  equation  (8): 

RDR)  +  RDitf  +  RDR f  =  RDR*A  +  A^RDR?  +  2k/, 

R^RD  +  DR^R  +  D  =  DR^  AR  +  R*  A^  RD  +  2k/ .  (10) 

Since 

f  (R^R)  =  &R  +  R)R  =  0 
at 

both  R^  R  and  R  '  R  are  antisymmetric  matrices.  Thus  in  equation  (10)  the  first  two  terms 
on  the  left  do  not  have  any  diagonal  elements.  Denoting  A  =  R^  AR  as  the  rotated  velocity 
gradient  matrix  we  have  the  following  equation  for  the  diagonal  elements, 

B%i  =  2A uDa  +  2k  , 

where  repeated  indices  do  not  indicate  summation.  If  we  assume  Da  =  e2pi  (pi  >  0  denotes 
extension  whereas  pi  <  0  represents  compression)  then  one  has 

Pi  =  An  +  ne~2pi  .  (11) 

Without  loss  of  generality,  we  arrange  the  eigenvalues  e2pl ,  e2p 2 , . . . ,  e2pd  in  non-increasing 
order: 

Pi  >  P2  ■  ■  ■  >  Pd  (12) 

where  d  is  the  dimension  of  space. 
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Extension 
e.  2"  «1 


Compression 

e2p<»\ 


Figure  2:  The  diffusive  term  (ne  2pi)  in  the  equation  for  pi  would  be  negligible  along  the 
extensional  axis  compared  to  the  compressional  direction. 


Having  worked  with  the  diagonal  matrix,  D ,  let  us  focus  on  the  eigenvector  matrix,  R. 
Considering  the  off-diagonal  terms  of  equation  (10), 

[Rj  RD]ij  =  [R^  R]n>D£j  =  [R^  R]ijDjj 
[DR^R]ij  =  Dn[R^R]£j  =  ~[R)  R\ijDa 

where  repeated  indices  do  not  indicate  summation.  Let  us  define  RJ  R  =  Q,  an  antisym¬ 
metric  matrix.  Thus  we  have, 


( Rjj  —  Dii-Aij  T  AjiDjj 

e2piAij  +  e2pi  Aji 


H ij  — 


pj  —  g2  pi 


i  +  j 


(13) 


Let  us  assume  there  exists  scale  separation  in  p’s  (pi  P2  Pd)-1  Thus  if  pt  pj, 

the  t  1  solution  is 


n 


*? 


_ 4. . 
A  ■■ 

3%i 


i  <  j 
i  >  j 


(14) 


Thus  A  is  independent  of  the  eigenvalues  and  we  can  solve  equation  (11)  directly, 


Pi(t)  =  pi( 0)  +  Ai(t)  +  \  log 


where  Ai(t)  =  /  Ah(t)cIt  . 


When  diffusion  vanishes  (k  =  0),  this  reduces  to 


1  +  2ne  2pi(0)  /  e  2Aiif') dt ' 

Jo 


(15) 


Pi(t)  =  pi( 0)  +  /  Ah(t)(It  . 

Jo 


(16) 


1This  is  almost  always  true  for  long  times,  but  it  may  fail  locally  or  due  to  a  symmetry  in  the  problem. 
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The  k  =  0  case  can  be  used  to  compute  Lyapunov  exponents2  from  the  eigenvalues  of  M: 


A  i  =  lim 


log(ePi^) 

t 


lim 

£—>•00 


Pi(t) 

t 


(17) 


The  convergence  of  the  above  limit  is  assured  by  Oseledec’s  multiplicative  ergodic  theo¬ 
rem  [7]. 

A  blob  of  initial  size  L  will  become  an  ellipsoid  with  the  length  of  its  main  axis  changing 
as  exp(Aji)  along  directions  corresponding  to  positive  and  negative  Lyapunov  exponents. 
The  rate  of  decrease  of  the  smallest  dimension  is  determined  by  A^-  Because  of  diffu¬ 
sion,  the  contraction  terminates  after  a  time,  tK  «  | A^| — 1  log(772 1 A^ | / ac) ,  when  the  diffusive 
lengthscale,  rK  =  \J K/\Xd\  is  attained. 


4  Large-deviation  form  of  the  central  limit  theorem3 

To  form  a  probabilistic  theory  of  mixing,  we  make  A,  and  therefore  pi  and  p-2  (we  restrict 
to  d  =  2  dimensions  from  now),  into  random  variables  and  deduce  the  approximate  long¬ 
time  moments  of  concentration.  To  do  so,  we  require  a  central  limit  theorem  (CLT)  which 
is  also  valid  for  large  deviations. 

Let  Xi  be  i.i.d.  (independent  and  identically  distributed)  random  variables  with  finite 
mean  x  and  variance  a2.  Let  x  be  the  mean  of  n  such  x.L.  Roughly  speaking,  the  classical 
CLT  asserts  that 

\  /  n(x  —  x')2  A 

P(x,  n)  ~  .  exp - ^ -  when  n  3>  1  and  x  —  x  <  -7=. 

v  '  V^W2  V  2  a2  )  ^ 

The  theorem  applies  to  large  n,  but  the  condition  that  x—x  <  ^  means  that  as  n  gets  large, 
the  region  around  the  mean  to  which  the  theorem  applies  grows  increasingly  narrow.  If  one 
is  concerned  about  the  tails  of  the  distribution,  the  theorem  is  not  powerful  enough.  This 
is  not  surprising  since  it  contains  only  second-moment  information  about  the  Xj  variables, 
and  higher  moments  may  be  critical  at  the  tails.  In  the  theory  of  mixing,  the  tails  of  the 
p  PDF’s  are  indeed  significant,  so  we  must  use  the  large  deviation  form  of  the  CLT,  which 
is  stronger  and  uses  information  from  all  moments  of  the  Xi  variables.  The  significance  of 
the  tails  will  be  evident  a  posteriori  from  the  analysis  using  the  large  deviation  CLT.  The 
large  deviation  CLT  asserts  that 

P(x,  n)  ~  exp  (— nS(x  —  x)) , 

where  S  is  the  Cramer  function  (or  rate  function,  or  entropy  function)  determined  by  the 
PDF  p(xi).  (In  these  notes  we  typically  neglect  algebraic  prefactors  in  asymptotic  relations, 
such  as  the  one  above.)  Whatever  the  form  of  p(xi),  S  is  convex,  and  5(0)  and  5'(0)  vanish. 
Taylor  expanding  5  about  0  in  the  large  deviation  CLT,  the  first  nonzero  term  is  quadratic 

2Lyapunov  exponents  describe  the  rate  at  which  infinitesimally  close  trajectories  in  a  dynamical  system 
diverge  from  each  other.  A  chaotic  system  is  usually  indicated  by  a  positive  infinite-time  Lyapunov  exponent. 

3The  ‘facts’  invoked  in  this  section  come  from  theorems  which  for  the  most  part  give  bounds,  not  equali¬ 
ties.  We  are  leaving  out  many  details  that  are  necessary  to  make  the  theory  rigorous.  See  for  example  [3,8,11] 
for  more  details. 
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in  x  —  x.  Truncating  after  this  term  recovers  the  classical  form  of  the  CLT,  so  the  classical 
form  simply  approximates  the  large  deviation  form  to  the  first  nontrivial  order.  The  explicit 
expression  for  the  Cramer  function  in  terms  of  p(xi)  involves  forward  and  inverse  Fourier 
transforms,  which  generally  cannot  be  explicitly  evaluated.  Rather  than  prove  the  large- 
deviation  CLT  or  derive  the  general  expression  for  the  Cramer  function  (see  e.g.  [3,8,10,11]), 
we  shall  merely  derive  the  Cramer  function  explicitly  in  the  simple  case  where  the  Xi  have 
a  Bernoulli  distribution,  though  the  general  formula  for  S  may  be  derived  by  an  analogous 
procedure. 


5  Large  deviation  CLT  for  the  Bernoulli  distribution 

Let  Xi  have  the  PDF  of  a  fair  coin, 

p(xi )  =  \6(xi  +  1)  +  \5(xi  -  1). 

The  characteristic  function  of  a  random  variable  is  the  Fourier  transform  of  its  PDF,  gen¬ 
erally  written  in  the  form  e~s^k\  We  wish  to  derive  p(x,n),  the  PDF  of  the  mean  x  of  n 
random  variables  a We  are  given  p(xi),  not  p(x,  n),  but  we  can  express  the  latter  in  terms 
of  the  former.  Because  the  xt  are  independent,  the  joint  distribution  of  x±,  ...,xn  is  simply 
the  product  of  their  individual  PDF’s,  so 

/n 

p{x  1,  ...,xn)s(x  -  \  y,  Xi)  dx  1  •  •  •  dxn 

i= l 

/n 

p(x i)  •  •  •  p(xn)5 (x  -  l  -  x)dxi  ■  ■  ■  dxn. 

i= i 

It  is  hard  to  compute  the  PDF  directly  from  the  PDF  of  Xi  by  the  above  formula.  Rather, 
we  shall  compute  the  generating  function  of  x,  and  then  its  PDF.  This  is  easier  because 
integrating  out  x-dependence  allows  us  to  then  integrate  out  the  Xi  variables.  The  generating 
function  of  x  is 

e-5(A.)  _  f  p^x^n^e-lkx(lx 


=  /  P(x i)  •  •  -p(xn)8 (x  -  \ 


Xj  e 


—ikx 


dx i • •  •  dxndx 


i— 1 


=  J p{x l)  •  •  •  p(xn)  exp  (  -  ^  x]dxi  •  •  •  dxr 

n  „ 

=  Y[  p(a 


i= 1 


ikx ^ 

[Xi  dxj 


1=1  ' 


p(x i)e  ™  dx i 

(s  J  +  1)  +  ""  l))e__srL dxi'j 


=  cos 
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We  recover  the  PDF  of  x  by  taking  the  inverse  Fourier  transform  of  its  generating  function: 

P(x,n)  =  ±  j  cos n(±)eikxdk 

=  i/exp(„ln[cos©]+to)* 

=  ^J  exp  (n[ln(cos  K)  +  iKxfjdK, 

where  K  =  K  The  last  form  of  the  above  integral  suggests  the  stationary  phase  approx¬ 
imation,  since  the  frequency  of  the  oscillatory  part  of  the  integrand  goes  to  infinity  as  n 
does.  Omitting  the  details,  the  stationary  phase  approximation  yields 

S(x)  ~  ~l(x  +  1)  In  (^f)  +  ln(l  -  x). 

The  classical  form  the  the  CLT  would  give  only 

p(x,n )  ~  exp  (  —  yS"(0)(x  —  x)2) 

~  eXP  {  -  ~2JW~ )  ■ 

In  this  example,  the  classical  CLT  overestimates  the  weight  of  the  tails  of  p(x,n),  though 
typically  the  opposite  is  the  case.  In  any  case,  the  large  deviation  form  is  needed  for  accuracy 
at  the  tails. 

6  Large  deviation  CLT  applied  to  mixing 

Recall  that  p2  has  generally  become  quite  small  at  long  times.  However,  in  a  fluid  with 
diffusion,  the  p-2  dimension  eventually  becomes  small  enough  that  further  contraction  by 
stretching  is  balanced  by  diffusion.  This  is  the  Batchelor  width,  \J ac/Ai,  beyond  which  the 
scalar  blob  will  contract  no  further.  As  such  late  times  we  can  assume  that  the  PDF  of  p2 
is  stationary  in  time,  so  as  t  — >  oo, 

P(pi,P2,t)  ~  exp  (  -  ts(^  -  Ai))pstat(P2)- 

The  exact  stationary  distribution  of  p2  is  irrelevant  to  the  heuristics  of  this  theory.  Suppose 
that  the  concentration  at  each  point  in  a  blob  of  scalar  is  inversely  proportional  to  the 
blob’s  total  volume,  so 

0(x,t)  ~  (det  M)1/2  =  e"(pl+p2). 

This  is  an  accurate  assumption  for  a  Gaussian  blob,  for  instance.  We  now  combine  the 
above  approximations  to  estimate  the  expectation  of  9a . 

(0a)(t)  ~  C(t)  j  e-“(P!+P2)  exp  ^  _  tS(^  -  Ai) Pstat(p2)dpidp2 

~  J  e~api  exp  —  —  X\)^jdpi, 
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Figure  3:  Generic  S(x)  and  S(h  —  Ai).  When  a  >  ac,  there  is  no  positive  h*  such  that 
S'(h*  —  A)  =  — a ,  so  we  must  choose  h*  =  0  for  the  saddle-point  approximation. 

where  C (t)  is  a  normalization  that  is  algebraic  in  t  and  is  therefore  asymptotically  dominated 
by  the  integral,  so  we  have  ignored  it  along  with  whatever  constant  emerges  from  integrating 
over  p2-  Letting  h  =  ^L, 

(6°){t)  ~  J  e~t(<ah+S{h-Xl))  dh  =  J  e~tHWdh , 

where  H(h)  =  ah  +  S(h  —  A).  This  integral  may  be  approximated  by  the  saddle-point 
method  for  large  t. 

The  saddle-point  method  typically  asserts  that 

(9a){t)  ~  e-W*, 

where  h*  is  such  that  H'(h*)  =  0,  i.e.  such  that  S'(h*  —  Ai)  =  —a.  In  the  present  case,  it  is 
not  always  possible  to  choose  such  an  h*  because  h  must  be  nonnegative.  This  is  clear  from 
the  generic  S  and  S'(h  —  A)  curves  plotted  in  Figure  6.  Whatever  the  exact  form  of  the 
Cramer  function,  S' (h)  begins  at  the  origin  and  increases  monotonically  in  h,  so  S' (h  —  A) 
begins  at  some  negative  point  (0,  —  ac)  and  increases  monotonically  in  h.  When  a  <  ac, 
H'(h*)  =  0  is  possible,  but  when  a  >  ac,  the  best  we  can  do  is  to  choose  h*  =  0.  That  is, 

(9a){t)  ~  e-t{ah*+s'(h*~x^) , 

where  h*  =  max{0,  A  +  (<S")_1(— a)}.  Neglecting  large  deviations,  the  classical  CLT  would 
have  given  a  decay  rate  of  — aA,  which  is  an  upper  bound  on  the  decay  rate  we  derived.  In 
fact,  the  tails  of  the  PDF  of  pi  act  significantly  to  slow  the  decay  of  concentration.  Above 
ac ,  ( 9a)(t )  is  constant  in  a.  We  may  understand  this  plateauing  as  high  concentration 
blobs  dominating  the  expectation  for  large  a.  These  are  the  blobs  that  have  experienced  no 
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stretching,  so  a  >  ac  moments  are  determined  solely  by  realizations  of  6  without  stretching. 
Note  that  the  plateauing  of  moments  also  means  that  the  La  norm  of  6  vanishes  as  a  — >  oo, 
which  is  characteristic  of  intermittency. 
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